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SUMMARY 

This paper documents the derivation of the data reduction equations for the calibration of 
the six-component thrust stand located in the CE-22 Advanced Nozzle Test Facility. The 
purpose of the calibration is to determine the first-order interactions between the axial, 
lateral, and vertical load cells (second-order interactions are assumed to be negligible). 

In an ideal system, the measurements made by the thrust stand along the three coordinate 
axes should be independent. For example, when a test article applies an axial force on the 
thrust stand, the axial load cells should measure the full magnitude of the force, while the 
off-axis load cells (lateral and vertical) should read zero. Fikewise, if a lateral force is 
applied, the lateral load cells should measure the entire force, while the axial and vertical 
load cells should read zero. However, in real-world systems, there may be interactions 
between the load cells. Through proper design of the thrust stand, these interactions can 
be minimized, but are hard to eliminate entirely. Therefore, the purpose of the thrust 
stand calibration is to account for these interactions, so that necessary corrections can be 
made during testing. These corrections can be expressed in the form of an interaction 
matrix, and this paper shows the derivation of the equations used to obtain the 
coefficients in this matrix. 


SYMBOLS 

CX 1 Aft-left axial calibration load cell 

C¥2 Aft-right axial calibration load cell 
CY1 Forward -right lateral calibration load cell 

CF2 Aft-right lateral calibration load cell 
CZ1 Forward-left vertical calibration load cell 
CZ2 Forward -right vertical calibration load cell 

CZ3 Aft-left vertical calibration load cell 
CZ4 Aft-right vertical calibration load cell 
RX 1 Aft-left axial reaction load cell 

RX1 Aft-right axial reaction load cell 
RY 1 Forward-left lateral reaction load cell 

RY2 Aft-left lateral reaction load cell 
RZ 1 Forward-left vertical reaction load cell 
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RZ2 Forward -right vertical reaction load cell 
RZi Aft-left vertical reaction load cell 
RZ4 Aft-right vertical reaction load cell 
FX Axial force, or force in the x-direction 
FY Lateral force, or force in the y-direction 
FZ Vertical force, or force in the z-direction 
MX Rolling moment, or moment about the x-axis 
MY Pitching moment, or moment about the y-axis 
MZ Yawing moment, or moment about the z-axis 

[Fmx\ The 5x1 force matrix that contains the rolling moment component 

[Fmy] The 5x1 force matrix that contains the pitching moment component 

[R] The 5x1 reaction load cell matrix that contains the rolling moment 

component 

[(?] The 5x1 reaction load cell matrix that contains the pitching moment 
component 

[*S] The 5x5 interaction matrix that contains the coefficients for the rolling 

moment calculation 

[L] The 5x5 interaction matrix that contains the coefficients for the pitching 
moment calculation 

Sij The coefficients in the matrix [5] 

Uij The coefficients in the matrix [ U] 

Gy The coefficients in the matrix [S] -1 

ju- The coefficients in the matrix [L] 1 

L P Forward moment arm for the pitching moment (distance from centroid to 

the CZ1 and CZ2 load cells) 

L P Aft moment arm for the pitching moment (distance from centroid to the 

CZ3 and CZ4 load cells) 

Ly Moment arm for the yawing moment (distance from centroid to the CY\ 

and CY2 load cells) 

Lr Moment arm for the rolling moment (distance from centroid to the CZ1 

and CZ3 load cells, and to the CZ2 and CZ4 load cells) 


THE CE-22 NOZZLE TEST FACILITY 

The CE-22 Advanced Nozzle Test Facility is located in the Engine Research Building 
(ERB) at the NASA Glenn Research Center. A schematic of the CE-22 Test Facility is 
shown in figure 1. The test section consists of a 23-ft long by T-Wfoot-inside-diameter 
vacuum tank. The upstream section is fixed while the downstream section slides open on 
rails, allowing for easy access to the test model. An inflatable rubber seal is used to 
prevent leakage when the two halves are closed for testing. 

The inlet air and exhaust is supplied by compressors and exhausters operated by Central 
Process Systems (CPS). The primary air system can supply inlet air at a pressure of 
40 psig with a maximum flow rate of 40 lbm/sec. There are also provisions for two one- 
inch secondary supply lines that enter through the top of the test tank. The secondary air 
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TO ALTITUDE/ 



Figure 1 . — The NASA Glenn Research Center Nozzle Test Facility (CE-22). 


can be supplied at 40, 150, or 450 psig with maximum flow rates of 21, 2, and 
10 lbm/sec, respectively. The exhaust system can pull a vacuum down to approximately 
1.9 psia to simulate altitude conditions at 48,000 feet. 

The CE-22 Test Facility also has a six-component thrust stand, which allows for 
simultaneous force and moment measurements along the three coordinate axes. The 
thrust stand and the altitude capability make the CE-22 Test Facility a unique asset to 
NASA. The facility is ideally suited for testing sub-scale advanced aircraft nozzle 
concepts which employ thrust vectoring. For further details about the CE-22 Test 
Facility, refer to references (1) and (2). 


THE SIX-COMPONENT THRUST STAND 

The installation of the thrust stand is shown in figure 2. The thrust stand is comprised of 
two major parts: the ground frame and the live frame. The ground frame is bolted to the 
floor of the test tank and the live frame is in turn attached to the ground frame through the 
8 reaction load cells. The inlet piping and the experimental model are mounted to the live 
frame. Thus, forces produced by the experimental model cause a displacement of the live 
frame relative to the ground frame. This displacement, on the order of thousandths of an 
inch, is converted to a voltage by each of the 8 reaction load cells. This voltage is then 
converted to a force reading by the data acquisition system. The interaction matrix is then 
applied by the data acquisition program to correct the individual force readings to take 
into account the interactions between the load cells. 
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Figure 2. — The six component thrust stand. 


There are three other items concerning the operation of the thrust stand that need 
mentioning: (1) There is a metric break that is located in the inlet air supply line just 
upstream of the thrust stand. The metric break is a physical break in the piping that 
isolates the experimental model and the part of the piping attached to the live frame of the 
thrust stand from the rest of the upstream air supply line. This break in the air supply line 
is enclosed by a labyrinth seal to minimize leakage; (2) The live frame of the thrust stand 
is actually comprised of two parts joined together by an elastic hinge. The purpose of this 
elastic hinge is to minimize the interaction between the axial and the lateral components. 
Because of the elastic hinge, off-center axial forces will not produce a yawing moment. A 
yawing moment can only be produced by offset lateral forces; and (3) the weight of the 
live frame, the piping, and the experimental model are zeroed out before the start of each 
test run. 

Figure 3 shows a schematic of the thrust stand. The sign convention used is as follows. 
The axial direction, or x-axis, is positive in the upstream direction. The lateral direction, 
or y-axis, is positive on the right hand side, aft-looking-forward. And the vertical 
direction, or z-axis, is positive in the downward direction. The terms used for the forces 
and moments are summarized in Table I. 
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Table I: Symbols and Terms Used for Forces and Moments. 


Symbol 

Force or Moment 

FX 

Axial force, or force in the x-direction 

FY 

Lateral force, or force in the y-direction 

FZ 

Vertical force, or force in the z-direction 

MX 

Rolling moment, or moment about the x-axis 

MY 

Pitching moment, or moment about the y-axis 

MZ 

Yawing moment, or moment about the z-axis 


LOAD CELL DESCRIPTION 

As previously mentioned, there are 8 reaction load cells that measure forces by detecting 
the displacement of the live frame relative to the ground frame. In addition, there are 8 
calibration load cells mounted in-line with the 8 reaction load cells. These calibration 
load cells are used to measure the applied load on the thrust stand during calibration. 

The calibration load cells are designated by a C (e.g., CXI, C71, etc.) and the reaction 
load cells with an R (e.g., RX 1, RY \ , etc). The positions of the 8 reaction load cells are 
summarized in Table II, and their exact locations relative to the centroid are shown in 
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Table II: Load Cell Locations. 


Load Cell 

Location 

RX\ 

Aft-Left 

RX 2 

Aft-Right 

RYl 

Forward-Left 

RY2 

Aft-Left 

RZl 

Forward-Left 

RZ2 

Forward-Right 

RZl 

Aft-Left 

RZ4 

Aft-Right 



figure 4. Forward and aft is relative to the direction of the airflow, and left and right is 
referenced from an aft-looking-forward perspective. 

Before the thrust stand itself is calibrated, each of the 16 load cells (8 calibration and 8 
reaction load cells) are sent out for calibration at the calibration laboratory. The load cells 
are calibrated against a standard that can be traced to a national standard. Each individual 
load cell works by converting a displacement into an electrical voltage. The calibration 
laboratory applies a known force on the load cell and plots a Force versus Voltage curve. 
This curve should be linear within the working range of the load cell. Once the slope and 
intercept from this plot are determined, these constants are entered into the data 
acquisition program. During an experiment, the displacement of the thrust stand is 
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measured as a voltage from the load cells, which in turn is converted to a force by the 
data acquisition system. 

A close-up of one of the vertical load cell pairs mounted on the thrust stand is shown in 
figure 5. The metal fixtures on each end of the load cell are flexures. Their job is to 
minimize the effect of off-axis forces on the load cell. Thus, ideally, an axial force 
produced by the experimental model is measured only through the axial load cells and 
should have little or no affect on the lateral and vertical load cells. 



Figure 5. — Calibration/reaction load cell pair. 


CALIBRATION PROCEDURE 

The purpose of the thrust stand calibration is to determine the first-order interactions 
between the load cells. For example, what effect does applying an axial force have on the 
vertical load cells? When a 100 lb. force in the x-direction is applied, one would expect 
the sum of the axial load cells to read 100 lbs. But due to constraints from the other 
attached load cells, this force may not be exactly 100 lbs. In addition, a force may be 
induced on the lateral and or vertical load cells. If the thrust stand is designed properly, 
such interactions can be minimized, but are hard to eliminate completely. The purpose of 
the calibration is to obtain a correction factor to account for these interactions. This 
correction factor is expressed as a matrix and is referred to interchangeably as the 
interaction matrix or the calibration matrix. 
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The calibration of the thrust stand is performed by applying known forces and moments 
and recording its effects on the reaction load cells. Forces are applied through the use of 
the stepper motors mounted in-line with the calibration load cells. The calibration load 
cells are used to determine the amount of force applied by the stepper motor. Forces and 
moments are applied in the manner shown in Table III. The first column indicates the 
type of load applied to the thrust stand. The second column indicates the calibration load 
cells used to measure this applied load. For the applied moments, the load is applied on 
one side of the thrust stand, then on the other side, as indicated by the "/". The third 
column shows the range and increment of the applied loads. The range and increment 
may be adjusted depending on the range of interest for upcoming test programs. 

As an example, when applying the moment about the y-axis, or MY, the stepper motors in 
line with CZ1 and CZ2 are activated to apply 0, -25, -50, -75, -100, -125, -100, -75, 
-50, -25, and 0 lbs. each to the thrust stand. The load is then applied in the other 
direction from 0 to 100 lbs. and back down to 0. At each point the readings indicated by 
the 8 reaction load cells are recorded. Readings are taken on the way up and on the way 
back down to determine the repeatability of the thrust stand. The same procedure is then 
repeated for the stepper motors mounted in-line with CZ3 and CZ4. 


Table III: Calibration Loadings. 


Applied Force or 
Moment 

Calibration Load Cell(s) used to 
Measure Applied Load 

Range and Increment 
of Applied Load 

FZ 

CXI, CX2 

0 to 3000, in. 250 lb 
increments 

FY 

CYl, CY2 

0 to -500 and 0 to 500, 
in 100 lb increments 

FZ 

CZ1, CZ2, CZ3, CZ4 

0 to -250 and 0 to 250, 
in 50 lb increments 

MX 

CZ1, CZ3 / CZ2, CZ4 

0 to -250 and 0 to 250, 
in 50 lb increments 

MY 

CZ1, CZ2 / CZ3, CZ4 

0 to -250 and 0 to 250, 
in 50 lb increments 

MZ 

CF1 / CF2 

0 to -500 and 0 to 500, 
in 100 lb increments 


The applied calibration loads are split equally between the two load cells when the FX, 
FY, MY, and MZ loads are applied. However, when the CZ1 and CZ3 (or CZ2 and CZ4) 
load cells are loaded together, as in the case when applying the FZ and MX loads, care 
must be taken in order to prevent an unwanted pitching moment. Due to the fact that the 
forward vertical load cells (CXI and CX2) have a different moment arm than the aft 
vertical load cells (CX3 and CX4), the forces must be applied such that: 

CZ1 = CZ2, (1) 


NASA/TM— 2003-212326 




CZ 3 = CZ 4 , 


( 2 ) 


and 


CZ 1 = 


L 


Pi 


L 


■CZ3, 


( 3 ) 


where L Pi = 22.25 in. is the distance along the x-axis from the centroid to CZ1 and CX2, 
and Lp 2 = 22.0 in. is the distance along the x-axis from the centroid to CZ3 and CZ4. 


THE FORCE AND MOMENT EQUATIONS 

The basic equation used to calculate forces and moments is: 

where [F] is a 6x1 matrix containing the 3 forces and 3 moments to be measured, [F] is 

the 8x1 matrix containing the 8 reaction load cell readings, and [5] is the interaction or 
calibration matrix. During the calibration, known forces and moments, [F], are applied to 
the thrust stand and the load cell readings, [F], are recorded from the data system. In 
order to obtain [5], its inverse, [S] -1 , must be determined from the following form of 
equation (4): 


[5]-'[F] = [F] (5) 

The problem with modeling the system using equations (4) and (5) is that [S] is a 6x8 
matrix, and is therefore not invertible. In other words, because there are 8 reaction load 
cells and only 6 components that need to be resolved, the system is indeterminate. The 
way around this problem is to divide the system into two parts. (This technique was 
developed by Roger Werner, 1991, NASA Glenn Research Center, OH, personal 
communication.) The load cells are combined, as shown below, so that there are 5 load 
cell readings and 5 components in each system. The two systems are represented below 
as: 


[iwHSH*] (6) 


~FX~ 



*12 

*13 

*14 

*15 

FX1 + FX2 

FY 


*21 

*22 

*23 

*24 

*25 

F71 

FZ 

= 

*31 

*32 

*33 

*34 

*35 

RY2 

MX 


*41 

*42 

*43 

*44 

*45 

FZ1 + FZ3 

MZ 


3i 

*52 

*53 

*54 

*55 _ 

FZ2+FZ4 
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[F my \ = [U][Q] 


'FX 1 


u n 

U 12 

U\ 3 

U u 

U 15 

RXI + RX2~ 

FY 


U 2 \ 

M 22 

w 23 

U 24 

U 25 

RYl 

FZ 

= 

M31 

W 32 

z/33 

W34 

M 35 

RY2 

MY 


W41 

W 42 

W43 

«44 

W 45 

RZl + RZl 

MZ 


_ M 51 

M 52 

M 53 

W 54 

M 55 _ 

RZ3 + RZ4 


( 8 ) 


( 9 ) 


In equation (7), the vertical load cells are combined left (RZl + RZl ) and right (RZl + 
RZ4) such that they produce a rolling moment, while in equation (9), the vertical load 
cells are combined forward (RZl + RZl) and aft (RZl + RZ4) such that they produce a 
pitching moment. See figure 3 for the load cell locations. 

In addition, the RX 1 and RX1 terms are combined, while the RY 1 and RY2 terms are kept 
separate. This is due to the elastic hinge. Because of this elastic hinge, an offset axial 
force would not produce a yawing moment; therefore, the terms can be combined. The 
RY 1 and RY2 terms need to remain separate in order to resolve the yawing moment, MZ. 

Both [Fmx\ and [F M y\ are usually calculated by the data acquisition program during an 
experiment. But because nozzles usually provide pitch and/or yaw vectoring, the 
components in \Fmy\ are more relevant since it contains the equation for the pitching 
moment. The calculation of [F M x\ serves as a check on [F MY \, and is also used to monitor 
the test article for unexpected rolling moments. 


THE INTERACTION MATRICES 

In order to better illustrate the function of the interaction matrix, the FX component as 
represented in equation (7) is written out as follows: 

FX = s n (RX 1 + RX 2) + s n RYl + s u RY2 + s l4 (RZl + RZl) + s l5 (RZ2 + RZ4) (10) 

In this form, it is easier to see what each of the coefficients represents. For instance, the 
coefficient sn represents the weighting factor for the effect of the sum of the two axial 
load cells, RX 1 and RX2, on FX. The coefficient s \2 represents the weighting factor for 
the effect of the RY 1 load cell on FX, and so on. Thus, for small interactions, one would 
expect the coefficient .s'u to be approximately one, and the other coefficients (in, s 13 , su, 
and si 5 ) to be close to zero. 

To get a better feel for the expected values of the coefficients in the interaction matrices 
[5] and [£/], assume for the moment that there are no interactions between the load cells. 
If this were the case, it is expected that the axial force would be the sum of the two axial 
load cells (iLTl + RX 2), the lateral force would be the sum of the two lateral load cells 
(i?71 + RY2), and the vertical force would be the sum of the four vertical load cells 
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(RZA + RZ2 + RZ3 + RZA). Similarly, the moments would be equal to the product of the 
applied force and its corresponding moment arm (see figure 4 for the lengths of the 
moment anns). For this ideal case, the following equations would be used to determine 
the three forces and three moments: 


FX=RXl + RX2 (11) 

FY = RYl + RY2 (12) 

FZ = RZl + RZ2 + RZ3 + RZ4 (13) 

MX = L R x (RZl + RZ3) + L R x (RZ2 + RZA) (14) 

MY = L Px x (RZl + RZ2) + L r , x (RZ3 + RZA) (15) 

MZ =L y xRY\ + L y xRY2 (16) 


These equations can be represented in the fonn of equations (7) and (9) as: 


' FX ] 


'1 

0 

0 

0 

0 " 

~RXl + RX2~ 

FY 


0 

1 

1 

0 

0 

RYl 

FZ 

= 

0 

0 

0 

1 

1 

RY2 

MX 


0 

0 

0 

L r 

Lr 

RZl + RZ3 

MZ 


0 

l y 

Ly 

0 

0 

RZ2 + RZA 


(17) 


FX~ 


'1 

0 

0 

0 

0 " 

~RX\ + RX2~ 

FY 


0 

1 

1 

0 

0 

RYl 

FZ 

= 

0 

0 

0 

1 

1 

RY2 

MY 


0 

0 

0 

L Pl 

L Pl 

RZI + RZ2 

MZ 


0 

Ly 

Ly 

0 

0 

RZ3 + RZA 


(18) 


Again, equations (17) and (18) represent the ideal case where there are no interactions 
between the load cells. However, for small interactions, the coefficients for the [5] and 
[f/] matrices obtained from the calibration should be fairly close to the above values. 
Thus the values in equations (17) and (18) can be used as a sanity check to determine if 
the calibration and data reduction were done correctly. 
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DERIVATION OF THE INTERACTION MATRIX COEFFICIENTS 


The previous sections explained how the interaction matrices, [5] and [U], are used to 
correct the forces and moments measured by the thrust stand. The following sections will 
explain the method by which the coefficient terms in the interaction matrices are derived. 

Equations (7) and (9), rewritten in the form of equation (5), are shown below: 


[S] 



(J \2 


<J \4 


" FX~ 


~RX\ + RX 2 ~ 

<J 21 

®22 


°24 

*^25 

FY 


RYl 

cr 3 i 

&32 

cr 33 

C34 

^35 

FZ 

- 

RY 2 

°41 

<J A 2 

°43 

C44 

°"45 

MX 


RZl + RZ 3 

_^Sl 

^52 

^53 

Cr 54 

^55. 

MZ 


RZ 2 + RZ 4 




m 

‘ \.F MY ] 

=m 


Mu 

M\ 2 

Mn 

M\ 4 

Ml 5 

FX 


~RX\ + RX 2“ 

Mn 

M22 

M23 

M2 a 

Ml 5 

FY 


RYl 


M32 

M33 

M34 

M35 

FZ 

= 

RY 2 

Mai 

Ma2 

MA 3 

M 44 

Mas 

MY 


RZl + RZ 2 

_Ms\ 

M52 

M53 

M54 

M55 _ 

MZ 


RZ 3 + RZ 4 


(19) 


( 20 ) 


( 21 ) 


( 22 ) 


Equations (20) and (22) are the basis for the calibration data reduction equations. The 
idea here is to apply a single known force or moment, for example, FX, record the 
reaction load cell readings, RX 1 + RX 1 , RY 1, etc., and determine the individual 
coefficients a and //. Once the individual coefficients are detennined, the [5] 1 and [t/]” 1 
matrices are inverted to obtain the interaction matrices [5] and [[/]. The next three 
sections explain in detail how this is done. Examples are given for the axial force (FX), 
the pitching moment (MY), and the yawing moment (MZ) calibrations. 1 The lateral force 
(FY) and vertical force (FZ) calibrations are similar to the axial force calibration, and the 
rolling moment (MX) calibration is similar to the pitching moment calibration. 


1 The term axial force calibration is used to indicate the loading of the axial calibration force. The axial 
force is not “calibrated” until all the interaction coefficients are determined from all the different loadings. 
Likewise for the terms lateral force calibration, pitching moment calibration, etc. 
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THE EQUATIONS FOR THE FX CALIBRATION 


For the axial force calibration, the stepper motors in line with the CXI and CXI load cells 
are used to apply the loads shown previously in Table III. For this case, all the terms in 
[Fmx\ are zero except for FX. Thus, equation (20) becomes: 



°"l2 

°"l3 

°"l4 

^15 

®2\ 

®22 

^23 

Cr 24 

°25 

°il 

^32 

C33 

C34 

^35 

o~ 41 

er 42 

O43 

cr 44 

°45 

&51 

°52 

^52 

°"54 

^55 


“ r fx~ 


~RXl + RX2~ 

0 


RYl 

0 

= 

RY2 

0 


RZl + RZ3 

_ 0 _ 


RZ2 + RZ4 


(23) 


From this, the first column of the inverse matrix can be detennined. Writing out the 
individual equations in the system, discarding the zero terms, and solving for era, 
equation (23) becomes: 


cr u FX = RX\ + RX2 

= > 

°ii 

RXI + RX2 
FX 

(24) 

<j 21 FX = RYl 

=> 

CT a 

RYl 
~ FX 

(25) 

<j 21 FX = RY2 

=> 

cr 3 i 

RY2 
~ FX 

(26) 

<j 41 FX = RZl + RZ2 

=> 

°41 

RZI + RZ3 
FX 

(27) 

a 5l FX = RZ2 + RZ4 

=> 

ff 51 

RZ2 + RZ4 
FX 

(28) 


Each of the era terms can be determined from the slope of the plot of the corresponding 
load cell readings versus FX. For example, the coefficient <j n is the value of the slope of 
the plot of RX 1 + RX2 versus FX. 


The coefficients jun, for the matrix [ U] 1 are found in a the same manner: 


Mn 


RXI + RX2 
FX 


Mu 


RY 1 
FX 


Mu 


RY 2 
FX 


(29) 

(30) 

(31) 
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RZI + RZ2 


(32) 


Max = 


FX 


Fsx 


RZ2+RZ4 

FX 


(33) 


Similarly, the FY and FZ calibrations will produce the coefficients <Ja and o# for [5] 
and the coefficients jUa and ju^ for [U] 


THE EQUATIONS FOR THE ME CALIBRATION 

The coefficients for the moment calibrations are determined in a similar fashion as the 
pure force calibrations, except there is an extra term due to the fact that an offset vertical 
force is used to produce the moment. For the ME calibration, the pitching moments are 
applied through two separate loadings. Forces are first applied through the stepper motors 
in line with the forward vertical load cells, CZ1 and CZ2. The procedure is then repeated 
with the aft vertical load cells, CZ3 and CZ4. The coefficients are computed separately 
for each loading, and then averaged in the final result. For both cases, the non-zero 
components of [F M y] are FZ and ME. Therefore, equation (22) becomes: 


Fu 

Fxi 

Fxi 

Fxa 

Fxs 

Fi\ 

Fll 

F23 

F24 

F25 

Fs 1 

Fil 

F33 

F34 

F35 

Fax 

Fa2 

Fa 3 

F44 

F45 

Fi\ 

F52 

Fs 3 

F54 

F55 


~~ 0 ' 


RX\ + RX2~ 

0 


RYl 

FZ 

- 

RY2 

MY 


RZl + RZ2 

_ 0 


RZ3 + RZ4 


(34) 


Writing out each of the individual equations, discarding the zero terms, and solving for 
fij 4 , equation (34) becomes: 


ju n FZ + <u u L p FZ = RX 1 + RX 2 


Fx4 


L 


1 ( RXI + RX2 3 
— Fn 


P v 


FZ 


/u 23 FZ + n 2A L p FZ = RY\ 
ju 33 FZ + /u 3A L p FZ =RY2 


/f>4 


L 


' p v 


RYl 

FZ 


■ // 23 


Fu 


L 


p v 


RY2 

~FZ~^j 


(35) 

(36) 

(37) 
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ju A3 FZ + n A4 L p FZ = RZ 1 + RZ 2 
jU 53 FZ + ju iA L p FZ = RZ3 + RZ4 



(38) 

(39) 


The jUj 3 terms in the above equations are determined from the FZ calibration and the 


RXI + RX2 RY 1 


etc.) are determined from the slope of the plot of the 


FZ FZ 

corresponding load cell readings versus the applied force, FZ. 


The 0/4 tenns are determined in a similar fashion through the MX calibration. 


THE EQUATIONS FOR THE MZ CALIBRATION 

For the MZ calibration, the yawing moments are applied through two separate loadings. 
Forces are first applied through the stepper motor in line with the CY 1 load cell. The 
procedure is then repeated for the CY2 load cell. The coefficients are computed separately 
for each loading, and then averaged in the final result. In this case, equation (20) 
becomes: 

a n cr, 2 cr 13 a 14 cr 15 0 RX1 + RX2 

rr 2 i ^22 ^23 ^24 ^25 FY RY] 

<t 3 i cr 32 ct 33 <t 34 cr 35 0 = RY 2 (40) 

(T 41 ^"42 ^ 43 ^"44 ^"45 ^ RZ 1 “h RZ 3 

0*51 °S2 ^53 ^54 [RZ2 + RZ4_ 

Writing out each of the individual equations, discarding the zero terms, and solving for 
O/ 5 , equation (40) becomes: 

a l2 FY + ct 15 L y FY = RX 1 + RX2 


u 22 FY + <j 25 L y FY = RY\ 


ct 32 FY + ct 35 L y FY = RY 2 
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^35 


L 


T V 


RY 2 
FY 


-<j 


32 


(43) 


a 42 FY + <j A 5 L y FY = RZl + RZ 3 


°~45 


L, 


1 f f?Zl + f?Z3 3 

— 1 <t 4 , 

Y V FT J 


(44) 


<j 52 FY + <j 55 L y FY = RZl + RZ 4 


1 f RZ2+RZ4 } 

55 L y { FY 52 J 


(45) 


The 0/2 terms in the above equations are determined from the FY calibration and the 

. RXI + RX2 RYl . , . , , . , , 

ratios ( , , etc.) are determined from the slope ol the plot ot the 

FY FY 

corresponding load cell readings verus the applied force, FY. 


The jUis terms for the [ U] 1 matrix are found in a similar fashion. 


SUMMARY OF THE INTERACTION COEFFICIENTS 

A summary of all the equations used for determining the coefficients in the [S] -1 and 
[f/] _1 matrices are shown in Tables IV and V. Once all the coefficients for the inverse 
matrices are found, the matrices can be inverted to obtain the interaction matrices [5] and 
[[/]. A sample set of calibration data is listed in Appendix A. 


Table IV: Coefficients cr for the [5] 1 Matrix 


i\j 

1 

2 

3 

4 

5 

1 

A(RX\+RX2) 

A(RX\+RX2) 

A(RXl+RX2) 

1 fA(RXl + RX2) ) 

1 { A(RXl+RX2) ) 

lyl AFY - ai2 J 

AFX 

A FY 

AFZ 

Lr l AFZ 13 J 

2 

A RYl 

ARYl 

ARYl 

1 f ARYl ) 

1 f ARYl ) 

AFX 

A FY 

AFZ 

Lr { AFZ a23 J 

Ly l AFY a22 j 

3 

ARY2 

ARY2 

ARY2 

1 f ARY2 ) 

1 f ARY2 ) 

AFX 

AFY 

AFZ 

Lr l AFZ 1733 J 

Ly l AFY 732 J 

4 

A(RZ\+RZ3) 

A(RZl+RZ3) 

A(RZl+RZ3) 

1 f A( RZ l+RZ 3) _ ) 

1 f At RZ l+RZ 3) ) 

rpvl af7 ^ 42 J 

AFX 

AFY 

AFZ 

Lr ( AFZ 43 J 

5 

A(RZ2+RZ4) 

A(RZ2+RZ4) 

A(RZ2+RZ4) 

1 (A(RZ2 + RZ4) ) 

1 f AIRZ2+RZ4) ) 

irl AF7 ~ CT52 J 

AFX 

AFY 

AFZ 

Lr l AFZ 53 J 
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Table V: Coefficients // (/ for the [£/] 1 Matrix 


i\j 

1 

2 

3 

4 

5 

1 

A(RX\+RX2) 

A(RX\+RX2) 

A(i«'l+«A'2) 

1 (A{RXURX2) 3 

MAUI-AU2) 3 


AFX 

A FY 

AFZ 

L P l AFZ m ) 

Ly l AFY nl ) 

9 

A RYl 

ARY\ 

ARY\ 

1 ( ARYl 3 

1 f ARYl 3 

z 

AFX 

A FY 

AFZ 

L P l AFZ M23 J 

Ly l AFY M22 ) 

9 

ARY2 

ARY2 

ARY2 

1 f ARY2 ) 

1 f ARY2 ) 

J 

AFX 

AFY 

AFZ 

Lp l AFZ /J33 > 

Ly l AFY M32 ) 

4 

A(RZ\+RZ2) 

A(RZl+RZ2) 

A(ffZl+«Z2) 

1 (A(RZURZ2) \ 

1 |A(«ZI-«Z2) ' 

i 

AFX 

AFY 

AFZ 

Lp i, AFZ 

Ly l AFY M42 j 

s 

A(RZ3+RZ4) 

A(RZ2+RZ4) 

A(i?Z3+«Z4) 

1 (A(RZ2 + RZ4) I 

l (A(RZ3 + RZ4) ) 


AFX 

AFY 

AFZ 

L P t AFZ ^ 53 J 

Ly l AFY M 52 ) 


CONCLUSION 

This report documents the derivation of the interaction matrix for the six-component 
thrust stand in the NASA Glenn Research Center’s CE-22 Advanced Nozzle Test 
Facility. It is intended to give the reader a better understanding of the theory behind the 
calibration and also insight into how the interaction coefficients are used during testing. 
Although the equations derived in the paper are specific to the thrust stand in the CE-22 
Test Facility, the theory and procedure can be applied to other similar thrust stands. 
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APPENDIX A— SAMPLE DATA 


Tables A1 to A9 show the raw data acquired from a typical calibration. The first column 
is the applied load, which is measured through the calibration load cells. The other 
columns are the reaction load cell readings for that particular applied load. All the load 
cell readings are in units of lbs-force. 

In order to reduce the data, each column of reaction load cell readings is plotted against 
the applied load. The slope of each plot is determined and is shown in bold below each 
column of reaction load cell readings. The slopes are then plugged into the appropriate 
equations in Tables IV and V to obtain the [5] 1 and [T/] 1 matrices as shown in 
equations (Al) and (A2), respectively. Finally, the matrices are inverted to obtain the 
interaction matrices, [5] and [U\, as shown in equations (A3) and (A4), respectively. 

NOTE: The slopes can be determined without plotting the data by using the SLOPE 
function from a spreadsheet application or by using some other curve fitting application. 
However, it is recommended that the plots be made to verify that the relationships are 
linear and that the load cell readings are repeatable. A non-linear relationship or non- 
repeating load cell readings may indicate possible impingement or other problems with 
the thrust stand. 
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Table Al. — Reaction Load Cell Data for an Applied Axial ( FX) Calibration Load 

FX RX1+RX2 RY1 RY2 RZ1+RZ3 RZ2+RZ4 RZ1+RZ2 RZ3+RZ4 

076 0.08116 -0.01616 0.09777 0.00806 -0.02446 -0.01633 -0.00007 
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Table A2. — Reaction Load Cell Data for an Applied Lateral ( FY) Calibration Load 

FY RX1+RX2 RY1 RY2 RZ1+RZ3 RZ2+RZ4 RZ1+RZ2 RZ3+RZ4 

173 1.72279 0.13895 -0.03260 0.16653 -0.56756 0.46330 -0.86432 
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Table A3. — Reaction Load Cell Data for an Applied Vertical ( FZ) Calibration Load 

FZ RX1+RX2 RY1 RY2 RZ1+RZ3 RZ2+RZ4 RZ1+RZ2 RZ3+RZ4 

182 0.21096 0.30270 0.60352 0.13023 -0.18184 0.03653 -0.08813988 
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Table A4. — Reaction Load Cell Data for an Applied Rolling ( MX) Calibration Moment (CZ2+CZ4 

FZ (CZ2+CZ4) RX1+RX2 RY1 RY2 RZ1+RZ3 RZ2+RZ4 

- 0.095 0.40605 - 0.01889 0.09239 0.07597 0.02983 
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Table A5. — Reaction Load Cell Data for an Applied Rolling ( MX) Calibration Moment (CZ1+CZ3 

FZ (CZ1+CZ3) RX1+RX2 RY1 RY2 RZ1+RZ3 RZ2+RZ4 

- 0.119 0.00983 0.04071 0.13592 - 0.02980 - 0.02729 
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Table A6. — Reaction Load Cell Data for an Applied Pitching ( MY) Calibration Moment ( CZ3+CZ4 

FZ (CZ3+CZ4) RX1+RX2 RY1 RY2 RZ1+RZ2 RZ3+RZ4 

0.114 0.08619 - 0.02446 - 0.06535 - 0.17562 - 0.02855 
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Table A7. — Reaction Load Cell Data for an Applied Pitching ( MY) Calibration Moment (CZ1+CZ2 

FZ (CZ1+CZ2 ) RX1+RX2 RY1 RY2 RZ1+RZ2 RZ3+RZ4 

- 0.057 - 0.08672 0.07834 0.13033 - 0.05664 - 0.02931 
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Slope 1 . 87379E-04 1.01435E-02 -2.36003E-03 9.96940E-01 9.95332E-04 



Table A8. — Reaction Load Cell Data for an Applied Yawing ( MZ) Calibration Moment (CF1) 

FY(CYl) RX1+RX2 RY1 RY2 RZ1+RZ3 RZ2+RZ4 RZ1+RZ2 RZ3+RZ4 

-0.060 0.34636 0.00540 0.09244 -0.14667 0.07093 -0.09759 0.02185 
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Table A9. — Reaction Load Cell Data for an Applied Yawing ( MZ) Calibration Moment (CY2) 

FY (CY2 ) RX1+RX2 RY1 RY2 RZ1+RZ3 RZ2+RZ4 RZ1+RZ2 RZ3+RZ4 

0.122 - 0.27797 - 0.00013 0.00000 - 0.19291 - 0.01109 - 0.14659 - 0.05741 
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